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O ■ ABSTRACT 

o 
o 

' Two-temperature spherical accretion flows produce ~ 100 Mev gamma-rays from the 

■ decay of neutral pions created in proton-proton collisions close to the black hole; they 
also produce ~ 10 keV X-rays by bremsstrahlung emission at large radii. The gamma- 

1^ , ray to X-ray luminosity ratio is nearly independent of black hole mass and accretion rate. 

CN i It does depend sensitively on the radial density profile of the accretion flow through the 



parameter a, where n ex r For the canonical Bondi value of a = 3/2, the gamma-ray 



J> ' to X-ray luminosity ratio is ~ 30. We interpret a recent Chandra detection coincident 



with the massive black hole at the Galactic Center as being thermal bremsstrahlung 



' emission from the accretion flow. With this normalization, the expected gamma-ray 



luminosity is ~ 10^^ ergs s ^ if a = 3/2. This is nearly two orders of magnitude 



m 
o , 

, above the detection threshold of the GLAST telescope. For a ~ 1/2, however, (a 

I value suggested by recent theoretical arguments), the expected gamma-ray luminosity 

' is only ~ 10^^ ergs s~^; GLAST should therefore provide an important probe of the 

I ' true accretion rate and radial density profile of the accretion fiow onto Sgr A*. 

■ Subject Headings: accretion, accretion disks — Galaxy: center — gamma rays: theory 



1. Introduction 

In roughly spherical accretion fiows, be it Bondi (1952) or advection-dominated accretion fiows 
(ADAFs; Rees et al. 1982, Narayan &; Yi 1994), the protons have temperatures comparable to their 
gravitational potential energy. Close to the black hole, a significant number of protons are energetic 
enough to exceed the threshold for the production of pions in proton-proton collisions. Neutral pions 
quickly decay to produce gamma-rays. It has long been recognized that this is a plausible source 
of gamma-ray emission from spherical accretion flows (Shvartsman 1971; Dahlbacka, Chapline, & 
Weaver 1974; Colpi, Maraschi, & Treves 1986; Mahadevan, Narayan, & Krolik 1997). 

In this paper we place the expected gamma-ray emission from spherical accretion flows on 
firmer observational ground by relating it to the more readily observable x-ray emission produced 
by thermal bremsstrahlung. A simple and relatively universal relationship between the two fiuxes 
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exists because both are produced by two-body processes. This is discussed in the next section (§2). 
We then apply these considerations to Sgr A*, the supermassive black hole at the center of our 
galaxy (§3). In §4 we briefly summarize our results. 

2. The Gamma-ray to X-ray Luminosity Ratio 

We first give a simple calculation of the gamma-ray to x-ray luminosity ratio assuming self- 
similar scalings for the density and temperature of the flow and a thermal distribution of protons. 
We then discuss the uncertainties introduced by these approximations. 

We take the temperature and number density of the flow to be 

Gp = ^0 and n = no r~°', (1) 

where 9p = kTp/mpC^ is the dimensionless proton temperature, r is the radius in the flow in units 
of the Schwarzschild radius (Rs), and no is the normalization of the density, which depends on, 
e.g., the black hole mass, the accretion rate, and the viscosity parameter a. For = 0.15, the 
proton temperature profile is that of non-relativistic Bondi accretion with an adiabatic index of 
7 = 5/3;0 comparable maximal temperatures and identical radial scalings occur in relativistic Bondi 
accretion (Shapiro 1973) and ADAFs (Narayan & Yi 1995). The electron temperature profile is 
rather uncertain; fortunately we will only need the electron temperature at large radii, r ^ 10^, 
where the flow is well approximated as one temperature. In equation (|l]) we allow the radial density 
profile to differ from the canonical Bondi value of a = 3/2; recent work on ADAFs has shown that 
much smaller values, e.g., a = 1/2, may be appropriate (see §2.2). 

The number of ~ 100 Mev gamma-rays produced per second and per cm^ is given by n^R{ep), 
where R{Op) is the reaction coefficient for thermal protons of temperature 9p. At 6p ~ 0.15, R{0p) 
can be approximated by R{9p) ~ i?o(^p/0.15)^, where Rq ^ 2x 10~^^ cm^ (see Fig. 3 of Dermer 
1986); for 9p ^ 0.05, R{9p) decreases much more rapidly than oc Op. Integrating over the flow, the 
photon luminosity in sa 100 Mev gamma-rays is 

N, = 4vri2|n^ R{ey-^- « ^i?|ngi?o. (2) 

Spherical accretion flows produce x-rays by thermal bremsstrahlung and by Comptonizing 
synchrotron photons. We are interested in very low luminosity systems where the former is expected 
to dominate. Thermal bremsstrahlung emission can also be expressed using a reaction coefficient. 
The number of x-rays of frequency v produced per second and per cm^ is given by n?R^(6e), where 
Rui9e) ~ pe7^^^exp[-hiy/kTe] and p ^ 1.3 x lO-^^ cm^ s"^ (e.g., Rybicki k Lightman 1979). 



^Relativistic corrections are small. 
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Integrating over the flow, the photon luminosity in x-rays of frequency v is 

Nx = AnRlnlpJ^ (^^^ r^-^^^e^^/^ eM-hi^/kTe]. (3) 

Equation (^) shows that at a frequency z/ the x-ray emission is dominated by the largest radius 
which satisfies kT^. ^ hu. This is because r^v?Te increases with increasing radius. We focus 
on x-ray emission at ~ 10 keV which is dominated by emission from r ^ 10^ — 10'*. At these 
radii the flow is quite accurately approximated as one-temperature so we can substitute Te = Tp 
{9e = mp9p/me) into equation (^) and perform the integral 

Nx « ^:^RsnlP%'^'rl '-'^, (4) 

where = 6^/6^, By = hu/nipC^, and /?' = [j{me/mp)^/'^ 3 x 10"^^ cm^ s"^. 

Combining equations (Q) and (^), the ratio of the gamma-ray luminosity at energy ~ 100 
MeV to the x-ray luminosity at energy Ex is given by 

^ (E^\ / 3.5-2a \ , 
Lx ^ \Ex) \ a ) - 

where the last approximation takes a = 3/2. 

Equation (|5|) shows that, for the self-similar analysis of this subsection, the gamma-ray to x-ray 
luminosity ratio of the flow depends only on the radial density profile. Since both pion decay and 
bremsstrahlung involve two-body processes the luminosity ratio from any spherical shell depends 
only on the local temperature(s). The radial density profile enters because pions are only produced 
in interesting numbers very close to the black hole while the x-ray luminosity primarily originates 
from rather large radii. For a = 3/2, equation @ predicts ~ 30Lx, an observationally interest- 
ing number (§3), while for a = 1/2, the predicted gamma-ray luminosity is certainly undetectable, 




2.1. Uncertainties in Bondi-like models (a = 3/2) 

Several models of quasi-spherical accretion (Bondi, ADAF) predict a nearly free-fall radial 
velocity and, as a consequence, a = 3/2. In this case, the primary uncertainty in the L^/Lx 
estimate of the previous section is the proton temperature: Op could be smaller than ~ 0.1 near 
the black hole. A priori this is quite worrying because of the strong temperature dependence of 
the pion reaction rate. 

We do not believe that this uncertainty poses a serious threat to the estimate of equation (|5|). 
Relativistic corrections to the classical Bondi solution are small (Shapiro 1973). As we explain 
below, the corrections due to rotation in the ADAF solution are relatively small as well. 
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In principle, ADAF models can have low proton temperatures near the event horizon (e.g., 
Op ^ 0.03). The low-temperature solutions, however, require small values of the dimensionless 
viscosity a, while numerical simulations and theoretical arguments (see §2.2) show that canonical 
ADAF models are only realizable if a is relatively large, roughly a ^ 0.1. 

For large a equation (|l|) is a reasonable approximation of even general relativistic calculations 
of the structure of ADAFs (Gammie & Popham 1998, Popham &; Gammie 1998; hereafter GP).^ 
For non-rotating black holes and a ^ 0.1, for example, our temperature profiles match those of 
GP very well;0they find maximal temperatures of Op ~ 0.1, consistent with our value. For rapidly 
spinning black holes, their temperatures are yet higher, reaching Op ~ 0.3. In addition, the density 
profile given by equation ([l|) is a reasonable approximation of the global calculations for large a. 
Self-similar solutions predict radial velocities ~ acs, where Cs is the sound speed of the gas. At 
small radii, however, the accreting gas must pass through a sonic point on its way into the black 
hole. For large a the "natural" radial velocity of the flow is of order the sound speed, so little 
deviation from self-similarity is required to match onto the sonic transition. 

To test the estimate of equation (|5|) we calculated the expected gamma-ray to x-ray luminosity 
ratio using several of GP's models and found generally good agreement. For a non-rotating black 
hole and an accretion flow with a ~ 0.3, for example, the more detailed calculation yields ~ 
lOLx for Ex = 10 keV, in reasonable agreement with equation (|5|). 

It is also important to emphasize that gamma-ray emission from pion decay is unlikely to be 
as sensitive to temperature as suggested by the simple thermal model we have considered. The 
collisionless plasmas of interest should efficiently accelerate protons to relativistic energies; in this 
case the total gamma-ray luminosity varies only linearly with changes in the thermal energy of the 
protons since there are always a substantial number of protons above the pion production threshold 
(see Mahadevan et al. 1997) .Q 

Equation (|5|) predicts a detectable gamma-ray flux only if a ?a 3/2. The above considerations 
suggest that equation (^) should be a good approximation in this limit. 

2.2. Non Bondi-like accretion models (a < 3/2) 

Modern theories and numerical simulations of quasi-spherical accretion flows suggest that the 
mean infall velocity can deviate substantially from the free-fall value, resulting in a < 3/2. There 

^This is not true for small a (e.g., Narayan, Kato, & Honma 1997). 

^GP consider several adiabatic indices for the flow; we compare only with 7 ~ 5/3, appropriate for a flow dominated 
by the energy density of the nearly non-relativistic protons. 

^Gruzinov & Quataert (1999) describe a proton heating model which yields very little proton acceleration. How- 
ever, if shocks or reconnection events occur in the accretion flow, a fraction of protons should be accelerated. 
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are three scenarios: convection-dominated accretion flows (CDAFs), winds, and turbulent heat 
conduction. 

CDAF: Two independent groups have performed numerical simulations of quasi-spherical non- 
radiating accretion flows with small values of the viscosity parameter (Stone, Pringle, & Begelman 
1999; Igumcnshchev & Abramowicz 1999). They both find a = 1/2 rather than a = 3/2. Narayan, 
Igumenshchev, &: Abramowicz (2000) and Quataert &: Gruzinov (2000) have explained this in terms 
of a CDAF. In such a flow angular momentum is efficiently transported inwards by strong radial 
convection. This nearly cancels the outward transport by magnetic fields, leading to a substantially 
suppressed accretion rate and a much flatter radial density profile. 

Winds: For large a ^ 0.1, CDAFs do not appear to be found in numerical simulations (Igu- 
menshchev Sz Abramowicz 1999); this is because the infall time of the gas is shorter than the 
convective turnover time, so convection is less dynamically important. For large a, however, a may 
differ from 3/2 for a different physical reason; strong outflows may drive away most of the accreting 
mass (Blandford & Begelman 1999; Igumcnshchev &; Abramowicz 1999). 

Turbulent heat conduction: Conduction preheats the infalling gas, reducing the accretion rate 
and flattening the density profile (Gruzinov 1999). 

3. Application to the Galactic Center 

Chandra observations of the Galactic Center detect a point source coincident with the non- 
thermal radio source Sgr A* to within 0.5" lO^Rs- Its 0.1-10 keV luminosity is Lx ^ 4:X 10^3 
ergs (Baganoff et al. 2000). It is very plausible that this represents the first x-ray detection of 
the supermassive black hole at the center of our galaxy. 

It is natural to interpret Chandra's detection as thermal bremsstrahlung from large radii in 
the accretion flow. As shown in §2, such emission would arise from r ~ 10^; the density required to 

match the observed luminosity is then ~ 4 x 10'^ cm^'^. The corresponding accretion rate is ~ 10~^ 
Mq yr~^, if the radial velocity of the gas is of order the sound speed. This is in good agreement with 
estimates based on the mass- losing stars in the central parsec of the Galactic Center (e.g., Coker 
&; Melia 1997; Quataert, Narayan, &: Reid 1999). The bremsstrahlung interpretation predicts the 
absence of short timescale variability in the observed x-rays (since the emission arises from large 
radii). It can also be tested by looking for x-ray line emission in deeper Chandra exposures (Narayan 
& Raymond 1999). 

EGRET observations of the Galactic Center region detect a source (2EG 1746-2852) with 
Lj « 10^^ ergs s~^ and a power law spectrum extending from sa 100 MeV to f« 10 GeV (Merck et 
al. 1996); it appears to be point like within the ~ 1° resolution of the instrument. Mahadevan et 
al. (1997) interpreted this emission as arising from an ADAF around the black hole at the Galactic 
Center. In their more comprehensive models of Sgr A*, however, Narayan et al. (1998) were unable 
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to produce gamma-ray emission at the required levels and satisfy other observational constraints. 
Moreover, the observed spectrum of 2EG 1746-2852 looks very similar to that expected from cosmic 
rays colliding with a dense cloud of molecular hydrogen. 

Our calculation in §2 predicts the expected gamma-ray luminosity from the accretion flow 
given the x-ray luminosity in thermal bremsstrahlung. We believe that the Chandra observations 
of the Galactic Center provide this thermal bremsstrahlung luminosity. With this normalization, 
equation (P) predicts ~ 10^^ ergs s~^ for a = 3/2. This is well above the ~ 2 x 10^^ ergs 
s~^ detection threshold of the forthcoming Gamma-ray Large Area Space Telescope (GLAST).0In 
addition, GLAST's angular resolution is expected to be significantly better than that of EGRET 
(for, among other things, the express purpose of identifying unidentified EGRET sources). GLAST 
will therefore likely have the capability of distinguishing a gamma-ray counterpart of Sgr A* (if it 
indeed exists) from 2EG 1746-2852. 

4. Discussion 

We have argued that the ratio of the ~ 100 MeV gamma-ray luminosity to the ~ 10 keV 
x-ray luminosity of a spherical accretion flow depends primarily on its radial density profile (where 
n (X r^"). In particular, for canonical spherical accretion flow models with a = 3/2, ~ 30-Lx; 
for a < 3/2, < Lx- 

Our analysis predicts the ~ 100 MeV gamma-ray luminosity expected from the accretion flow 
onto the supermassive black hole at the center of our galaxy (§3): ~ 10^^ ergs s~^ if a = 3/2 
while ~ lO'^^ ergs s~^ if a = 1/2. For a = 3/2, this estimate is nearly two orders of magnitude 
above the detection threshold of the GLAST telescope. We expect, however, that no gamma-rays 
will be observed coincident with the black hole, supporting theoretical suggestions (see §2.2) that 
the density profile in spherical accretion flows is significantly flatter than the canonical r~^/^ profile. 

We thank Charles Gammie and Bob Popham for the use of their global models of ADAFs, Fred 
Baganoff for information about Chandra observations of the Galactic Center, and John Bahcall 
for discussions and comments. EQ is supported by NASA through Chandra Fellowship PF9- 
10008, awarded by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical 
Observatory for NASA under contract NAS 8-39073. AG was supported by the W. M. Keck 
Foundation and NSF PHY-9513835. 
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